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ABSTRACT 

Context. The discovery of a 2M Q neutron star provided a robust constraint for the theory of exotic dense matter, questioning the 
existence of strange baryons in the interiors of neutron stars. With many theories failing to reproduce this observational result, several 
equations of state containing hyperons are consistent with it. 

Aims. We study global properties of stars using equations of state containing hyperons, and compare them to those without hyperons 
in order to find similarities, differences and limits that can be compared with the astrophysical observations. 

Methods. Rotating, axisymmetric and stationary stellar configurations in General Relativity are obtained, and their global parameters 
are studied. 

Results. Approximate formula? describing the behavior of the maximum and minimum stellar mass, compactness, surface redshifts 
and moments of inertia as functions of spin frequency are provided. We also study the thin disk accretion and compare the spin-up 
evolution of stars with different moments of inertia. 
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^ 1. Introduction 



O The discovery of a 2M Q neutron star (NS) (Demorest et al. 



C/3 
C3 



2010) provides a strong motivation for the dense-matter physics 
community to understand the role of exotic phases in neutron 
j stars' interiors. The appearance of new particles (in general, 
new phases of matter) typically results in the softening of mat- 
^ — I ter, since the contribution to pressure from energetic nucleons is 
replaced by that of slowly-moving, massive particles; the pres- 
sure increases less steep with the density in the equation of state 
(EOS), which in turn lowers the NS maximum allowable mass 
L^j (for a review concerning dense matter EOSs, see |Haensel et al.| 
" . |2007| l. On the other hand, the state of the dense matter in the NS 
cores is very different than the one known at the energies and 
densities reached by the terrestrial experiments; notably, stable 
strange matter may appear in the interiors of NSs at densities a 
few times the nuclear saturation density, strange baryons - hy- 
>• perons - being a very probable result. 



X 



We study a set of selected EOSs, representing current the- 
oretical approaches to the description of the dense matter, that 
are compliant with the robust constraint put forward by the ob- 
servers. If we assume that hyperons do exist in the interiors of 
massive NSs, then based on our present knowledge it seems nec- 
essary for such objects to be composed of sufficiently stiff matter 
(at lower densities) in order to sustain the softening that is intro- 
duced by the appearance of hyperons - that is, substantial stiff- 
ness of the EOS of nucleonic (non-strange) matter is a necessary 
condition for the existence of hyperons in the core of a massive 
NS. Consequently, stellar configurations that contain hyperons 
will differ from those without hyperons. We therefore seek sim- 
ilarities and differences between them, as well as limits that can 
be applied to astrophysical observations. 



This note is composed as follows: Sect. [2] contains the de- 
scription of EOSs used. Sect. |3]presents the results for the grav- 
itational mass, surface redshift, moment of inertia and accretion 
tracks. Sect.[4]contains conclusions and summary. 



2. Equations of state 



Following selec tion of 'standard' EOSs is used: APR EOS ( |Ak-| 
mal et alj l998 model A18 + 6v + UIX*) is a variational, non- 
relativistic many-body solution with relativistic corrections. DH 
EOS (Douchin & Haensel 2001) is constructed using a non- 
relativistic energy density functional based on the SLy4 effective 
nuclear interaction, designed to describe both crust and core in 
an unified way. BSK20 EOS (moderate stiffness EOS by Goriely 
|et al.|[20T0| is also based on the nuclear energy-density func- 
tional theory, using the generalized Skyrme forces fitted to ex- 
perimental nuclear data and reproducing properties of infinite 
nuclear matter from the calculations of many -body interactions. 

Realistic microscopic dense-matter theories predict maxi- 
mum masses of NSs with hyperons much below the observed 
2M Q (see e.g., |Burgio et al.||20lT) . This problematic fea- 
ture is usually cured by providing stronger repulsion between 
the hyperons and/or letting them appear only at very high densi- 
ties. Below we list a selection of EOSs with hyperons, derived 
from theories that successfully deal with the problem of a mas- 
sive NS: DS08 EOS ( |Dexheimer & Schramm|2008) l uses an ef- 
fective hadronic SU(3) chiral model including the baryon octet 
and fourth-order self-interaction terms of the co, p and <p vector 
mesons. GM1Z0 EOS ( Weissenb orn et aL||2012| ) is based on a 
relativistic mean field model which allows studying the depar- 
ture from a vector meson-hyperon couplings stemming from the 
SU(6) quark model to a more general SU(3) prescription; the 
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Table 1. Selected parameters of static configurations (mass expressed 
in M Q , radius R in km, the moment of inertia / in 10 45 g • cm 2 and the 
frequency f* = (l/2n) yjGM^/RiM^Ji in Hz; R lA and I 1A are taken 
for.V/- 1.4 M ). 



EOS 



DH 



R 



1.4 



'1.4 



APR 2.19 
BSK20 2AT 



2.05 



GM1Z0 2.29 



TMIC 2.05 



9.93 
10.17 



2.24 

T00~ 



11.34 
11.75 



1.31 
1.39 



9.99 



12.51 



2740.47 
2630.50 



2.27 11.73 1.37 2630.63 



BM165 2.03 10.68 2.24 13.59 1.74 2367.82 
DS08 2.05 12.02 2.56 13.91 1.81 1993.42 



12.00 3.09 13.89 1.84 2108.51 



2.72 14.51 1.91 1873.41 



3.1. Gravitational mass 

Fig. [T] shows how the mass-central baryon density M(«b) rela- 
tion changes with / for selected EOSs. A strong reduction of 
available mass and central density range for large / is charac- 
teristic to hyperonic EOSs; moreover, configurations near the 
mass-shedding limit (left sides of curves) may be more massive 
than those near the axisymmetric perturbation instability limit 
(right sides). Fig. [2] shows the range of available masses as a 
function of spin frequency /; we note that the mass of a hyper- 
onic EOS star is already quite constrained by the existence of 
716 Hz pulsar ( [Hessels et aL|2006] i. Detection off 800 Hz 
pulsar would mean, that such a stiff EOS cannot yield a star with 
M < 1.4 M Q (in the case of the TMIC EOS this minimum mass 
is approached already for 716 Hz!). The growth of the maxi- 



GM1 model with the ratio between the meson octet and the sin- 
glet coupling constant z - g&/gi = (all the baryon-meson cou- 
plings are equal) was used, which makes it the stiffest in our sam- 
ple. BM165 EOS ( |Bednarek et al.|201l) was also derived from 
a relativistic mean field model, with a non-linear Lagrangian 
that includes quatric terms in meson fields, and two additional 
hidden-strangeness mesons cr* and <p that couple to hyperons 
only. TMIC EOS (Gusakov et al. 20121 employs a prescrip- 
tion similar to the BM165 EOS for purely nucleonic matter, and 
also introduces scalar cr* and vector <f> mesons in the description 
of hyperon interaction, with an additional A - A hyperon attrac- 
tion. SU(6) symmetry breaking applied in the model amounts to 
2 = 0.2. 

3. Results 

Following subsections contain the results for constant spin fre- 
quency tracks (sequences of configurations) of rotating stars. 
The frequency range spans astrophysically-relevant range from 
/ = Hz (static configurations, see Table [T] for their represen- 
tative parameters) up to / = 1200 Hz (i.e., much above the fre- 
quency of 716 Hz of the most rapid pulsar to date, PSR J1748- 
2446ad of |Hessels et al.|2006| l. Rigidly-rotating, stationary and 
axisymmetric stellar configurations were obtained by means of 
the numerical library LORENEQnrotstar code, using the for- 
mulation of |Bonazzola et al7| (|1993]>, with the accuracy ch ecked 
by a 2D virial theorem (IBonazzola & Gourgoulhon 1994). The 
sequences are limited by following conditions: from the low cen- 
tral density end, they terminate at the so-called mass-shedding 
limit, when the orbital frequency of a test particle at the stars' 
equator equals to the stellar spin frequency /. The high central 
density limit is marked by the onset of the axisymmetric per- 
turbation instability, described by the condition dM/dA\j = 0, 
where M is the gravitational mass and A a suitable parametriza- 
tion of the sequence of configurations (e.g., the central density; 
see |Friedman et al.|1988| for details). We define the moment of 
inertia as / = where J is the total stellar angular momen- 
tum and O is the angular frequency, D, = 2nf (for the definitions 
of M and 7, see Bonazz ola et al.|1993 1. For non-rotating config- 
urations, / is calculated using the slow-rotation approximation 



(Hartl e|1967) . 

Parameters of the approximate formulae presented below 
were obtained by the x 2 fitting, yielding a typical accuracy of 
the order of one per cent. 
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Fig. 1. (Color online) mass-central baryon density M - n b rela- 
tions for selected spin frequencies /, for stars composed of the BSK20 
EOS (dashed red lines) and three hyperonic EOSs: DS08 (solid green), 
BM165 (solid black) and GM1Z0 EOS (dashed blue lines). From bot- 
tom to top for each EOS, / equals 0, 900, 1000, 1 100 and 1200 Hz. 
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Fig. 2. (Color online) Range of allowable gravitational masses as a 
function of spin frequency /. Vertical line marks the 716 Hz frequency 
(spin frequency of the most-rapidly rotating NS to date). 

mum mass M max , normalized by the maximum mass of the static 
configuration M,^ ax can be approximated by 



^max/M^ ax = 0.49(/// S ) 8/3 + l, 



(1) 



http : // www . lor ene . obspm . f r 



whereP = (1/2^) y GMf wdx /R (M„ ax ) 3 is the frequency scaling 
factoirlformed of static maximal parameters, gravitational mass 

2 / s equals to the orbital frequency of a test particle orbiting spherical 
mass Mj^ ax at a distance R(M^). 
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M. Bejger: Parameters of rotating neutron stars with and without hyperons (RN) 
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Fig. 3. (Color online) Maximal mass M max (left panel) and mini- 
mal mass M ra j n (right panel) of rotating NSs as a function of the spin 
frequency /. Values are scaled by the parameters of non-rotating con- 
figurations. 



M^ ax and corresponding radius 7?(M^ ax ). In the case of mini- 
mum mass M m ; n , one can obtain a similar (although somewhat 
cruder) approximation to the one for the M max : 



M min /A4 ax = 3.57 (ft?) + 0.1. 



(2) 



Above relation underestimates the M min /M^ ax for BM165 EOS 
stars for a given frequency, because the maximum-mass radius 
/?(M max ) is smaller than for other hyperonic EOSs, resulting in a 
larger/ 8 . 

3.2. Surface redshifts 

Surface redshift may provide important information about the 
spacetime in vicinity of the NS. In the case of a rotating star, 
one considers the redshift of photons coming from the pole, Zp, 
as well as two equatorial redshifts for photons, emitted tangen- 
tially in and opposite the direction of rotation (forward Zf and 
backward Zb redshift, respectively; for definitions, see e.g., |Gour^1 
goulhon 2010, Sect. 4.6). Left panel in Fig.[4]shows the relations 
for z for the minimum compactness (M//? eq ) m ; n , i.e., configura- 
tions at the mass-shedding limit, that bear some resemblance to 
those during the photospheric radius expansion burst (Stei ner et] 
aL||2010 and references therein); R eq denotes the equatorial ra- 
dius). These relations may be regarded as upper limits on the 
compactness parameter for a measured z. We approximate z p 
and Zb by a (M/R) 3/2 + b, where a = 9.62 and b = 0.023 for 
Zp, and a = 25.11 and b = 0.18 for Zb', Zt is approximated by 
-0.68 (M/Kf A . 

Right panel shows redshift functions for the maximum- 
compactness (M//? eq ) max configuration along the line corre- 
sponding to the axisymmetric perturbation instability limit. For 
spin frequencies under consideration, it may be also treated as 
an approximate lower limit on the stellar mass. We propose 
to parametrize these relations as follows: for £ = (Zb - Zt)IZp, 
M/M s max * 0.008£ 25 + 1; for f = 1 - z f /z p , M/M^ ax is approxi- 



mated by 0.Q6C 



2.5 
15 + 1. 



In the case of the third function, we chose 



£ = 1 - z f /zb, and M/M, s nax a 0.04^ b + 0.02^ + 1. The value of 
Zb - zt quantifies the redshift effect on the maximal line broad- 
ening coming from the equator of a rotating star. Fig.|5]presents 
the maximal compactness for such stars against the two equa- 
torial redshifts - the (M//? eq ) max stay remarkably constant with 
rotation, and could be used to distinguish the two sets of models. 
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Fig. 4. (Color online) Surface redshifts for rotating configurations. 
Left panel: Equatorial backward (top curve) and forward (bottom curve) 
redshifts z b and z f , and the polar redshift z v for configurations with the 
minimal mass, as functions of the compactness (M//? e q)min- Right panel: 
Functions of the three redshifts for rotating, maximally-compact config- 
urations. 
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Fig. 5. (Color online) Maximum compactness vs equatorial redshifts 
Zf and Zb for different rotation rates; for Zt, f increases to the left, for 
Zb to the right from the bigger marker, denoting the static configuration. 
Points correspond to frequencies from to 1200 Hz. 



3.3. Moment of inertia 

A potentially interesting NS parameter is the moment of iner- 
tia / (in principle measurable using the spin-orbit coupling, |Lat-| 
|timer& Schutz||2005[ l. Fig. [6] shows the behavior of the grav- 
itational mass Mas a function of the moment of inertia / for 
static configurations. The curves in the left panel may be ap- 
proximated by a straight line, I(M) = a\M + b\. For normal 
EOSs, a\ = 1.22 and b\ = -0.34, for hyperonic EOSs a\ = 1.51. 
On the right panel we plot the value of I/R b as a function of M; 
this value enters the estimation of the minimum magnetic dipole 

field at the pulsar surface, B > (l/R^) U2 (3c 3 PP/8tt 2 ) 1/2 , where 

P = In jf denotes the spin periocjj Fig. [7] presents how the 
available range of the moment of inertia changes with rotation. 
For the frequencies / > 800 Hz the moment of inertia tend to 
be already quite constrained. The minimum value of / for sub- 
millisecond rotation coincides with the high-density end of the 
track, while the mass-shedding limit is characterized by a large 
/ (because of a strong equatorial radius dependence, the effect is 



3 Despite the fact that / grows with rotation, the value of I/R 6 de- 
creases because of the high power of R in the denominator. The static 
results are therefore upper limits for the value of I/R 6 for a given M. 
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much more pronounced than for the gravitational mass, Fig.[2]i. 
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Fig. 6. (Color online) Left panel: moment of inertia / as a function 
of gravitational mass M for static solutions. Right panel: moment of 
inertia / divided by R 6 (as used in the estimation of the minimum surface 
magnetic field of a pulsar) for static stars. Vertical lines correspond to a 
fiducial model of / = 10 45 g • cm 2 and 10 km. 
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Fig. 7. (Color online) Top panels: The behavior of I(M) for rotating 
stars (from left to right / equals 0, 900, 1000, 1100 and 1200 Hz) for 
hyperonic EOSs and DH EOS for comparison. Bottom panels: Allow- 
able range of the moment of inertia for (left panel for purely nucleonic 
stars, right panel for hyperonic stars). Thin (thick) lines correspond to 
the minimum (maximum) of / that a stable configuration can attain at a 
given frequency. 

As many softer EOSs are now ruled out by observations, we re- 
produce the formula suggested byJLattimer & Pr akash (200T|) in 



the form of Bejger & Haensel (2002), with slightly different co- 
efficients in the static case. For an astrophysically relevant range 
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Fig. 8. (Color online) Left panel: relation between I /MR 2 and the 
compactness parameter M/R for static configurations (grey) and those 
rotating at 1200 Hz, in which case R = R cq . Right panel: Minimal 
M/R eq ratio (value corresponding to the mass-shedding configuration) 
as a function of rotational frequency. 



of masses, M > 0.5 M Q , we have 
I /MR 2 = fl (M/R) + b (h 



(3) 



with «o = 1-40 and bo = 0.19 (left panel in Fig. [8}. This approx- 
imation is generally valid for the static results in our sample. A 
more universal formula, which takes into account rotation, will 
depend on the spin frequency /; the compactness M/R must be 
re-defined as the gravitational mass M to the circumferential ra- 
dius R eq ratio. For EOSs without hyperons, we propose a follow- 
ing approximation: 

a(f) = 0.l9f*+a 



and b(f) = -0.042/f + bo- 



(4) 



where f\ - f/(l kHz). Since the rotation influences hyperonic 
stars in a different way, we produce a analogous approximation 
for a subset of hyperonic EOSs only - in the static case we obtain 
aoh =1.14 and boh — 0.22, and for rotating configurations 



a(f) = 0.51 f?+a 0h and b(f) = -0.099/j + b { 



'0h- 



(5) 



In order to estimate the moment of inertia, one needs to know the 
minimum value of the compactness parameter M/R eq for a given 
spin frequency (right panel on Fig.lHll. It can be described by 



(m/rJ . =a c tf + b e , 



(6) 



where b c = 0.005, and a c depends on the type of the EOS: 
for EOSs without hyperons a c = 0.074, otherwise a c = 0.099. 
Thanks to their stiffness, hyperonic EOSs produce larger radii 
near the mass-shedding limit, but the minimal mass for a given 
/ is also higher in comparison to standard EOSs (see Fig. [2J, 
which results in higher minimal M/R eq . 

3.4. Spin-up by disk accretion 

A process, in which the moment of inertia is relevant is the so- 
called recycling of pulsars to millisecond periods. We use the 
Bejger et al.|(j2011a[) - m agnetic torque in 



model described in 



the form given by Kluzniak & Rappaport 2007 improved for the 
existence of the marginally-stable orbit, and the magnetic field 
decay proportional to the amount of accreted mass - to study the 
relation of spin-up to the EOS. The evolution of the total stellar 
angular momentum in the process of accretion is 



dJ/dM b = 1-1 



mag 



(7) 
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M. Bejger: Parameters of rotating neutron stars with and without hyperons (RN) 



where / is the orbital angular momentum of a particle in the disk 
per unit baryon mass, and Z mag is the contribution from the brak- 
ing magnetic torque, resulting from the magnetosphere-disk in- 
teraction. 

For comparison with the realistic EOSs, we employ an use- 
ful fiducial modefll45: a star with a constant radius R = 10 
km, and a constant moment of inertia I = 10 45 g ■ cm 2 . To 
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Fig. 9. (Color online) Spin evolution of an accreting star as a function 
of the accreted mass (accretion rate is M = 10~ 9 M Q /yr, the horizontal 
axes correspond to accretion time f acc in Gyr) with initial mass M m [ = 



1.4 M Q . Left panel: accretion with the magnetic field decay, B n 
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Fig. 10. (Color online) As in the left panel of Fig. [5] but for M ini = 
1.2 M (left panel) and M ini = 1.9 M (right panel). 

compare the effect of the EOS on the recycling process, Fig. [9] 
shows the evolution of spin frequency / for a chosen accretion 
rate, M = 1(T 9 M Q /yr and initial magnetic field B = 10 12 G (left 
panel). As anticipated, stars with larger moments of inertia gain 
/ less rapidly: with the above parameters, a hyperonic star has 
to accrete almost twice as much matter as a 'standard' EOS star 
to reach given / (see also Be jger et al.|201 la) . For comparison, 
right panel shows the importance of the magnetic torque - when 



4 In general, the changes of stellar parameters relate to each other as 
dM = CldJ + /udM b , where ^ = l/u 1 is the chemical potential per unit 
mass, transfered onto the star by the accreting particle, u' being the 
time component of the particle 4-velocity. In the case of B = 0, n = 
e - m, where e is the specific energy of a particle infalling onto the 
star (Sibgatullin & Sunyaev 2000, Zdunik et al. 2002), the gravitational 
mass growth of such a star is therefore dM = edMt\ we adopt this 
prescription for B t too (we also correct a misprint in Eqs. 4 and 5 of 
Bejger et al. 201 la instead of u', it should read l/u'). 



it is neglected (B = 0), the amount of mass needed to spin up the 
star to a given / is much smaller (lines end at the mass-shedding 
limit). By rewriting Eq.[7]for J = IQ., 



dQ/dM ft = (/tot - QdZ/dMi,) //, 



(8) 



one deduces that both large I and its growth with M decrease the 
spin-up rate (to be compared with the fiducial 145 configuration; 
see also the top left panel in Fig. [6| where we compare BSK20 
and DS08 EOSs). We also note that, contrary to intuition, more 
massive stars are spinning up faster, as shown on Fig. [10] for 
= 1.2 M Q and Mj n ; = 1.9 M Q . This feature is an outcome 
of the fact that l tot is larger for larger M - in the general case 
the expression for k t(M) is rather complicated (see Bejger et 
[aL]|201 la] for details); in the simplest case, however, for B = 

in the Schwarzschild case, I oc -\jMr ms , and r ms - the radius 
of the innermost stable circular orbit - depends linearly on M. 
This effect is independent of the change of I, as shown for the 

1 = const, tracks with different M m \. The change of d//dM^ for 
realistic EOSs, related to the KM) behavior near the M max , is 
visible as a slight change of slope in the right panel of Fig. [10] 
just before the accretion ends at the instability limit. 



4. Summary and conclusions 

We have studied a set of recent EOSs that contain hyperons and 
yield static NSs consistent with a robust observational constraint, 
a 2M G NS mass measured by ( |Demorest et al.|2010| l. Stars con- 
structed using hyperonic EOSs were compared with these us- 
ing minimalistic approach i.e., containing nucleons only. Their 
global parameters: gravitational mass, compactness, surface red- 
shift and moments of inertia were calculated for a broad range 
of spin frequencies (0 - 1200 Hz). We have focused on extremal 
values of these parameters and find that the minimum mass for 
hyperonic EOSs increases much faster with rotation than in the 
case of standard EOSs - for frequencies just slightly larger than 
716 Hz the minimum mass approaches 1.4 M G . The same ap- 
plies to the compactness parameter, M/R ecj - its minimal value is 
substantially larger than for standard EOSs at a given / (Fig. [8}. 
We also note that the sub-millisecond rotation confines the mass 
into a narrow interval (a feature that is present for 'standard' EOS 
stars for much higher frequencies, see Bejger et al. 2007). 

A number of approximate formulae describing the whole set 
of EOSs is provided: we approximate the relations between the 
maximum and minimum mass, minimum compactness M/R eq 
and functions of surface redshift for minimal and maximal com- 
pactness configurations, for a range of studied frequencies. We 
find that the maximum compactness configuration has M/R eq 
that changes little with rotation, but separates the two groups. 

Due to their stiffness, hyperonic stars have comparatively 
large moments of inertia that increase further with rotation to 
values significantly larger than the usual 10 45 g • cm 2 - we extend 
the parametrization of Bejge r~& Haensel|200 2 by supplying an 
approximate 7/ (MR 2 ) vs M/R relation for rotating stars. Large 
moment of inertia is also one of the parameters in the process of 
"recycling" of a millisecond pulsar that hinders the spin-up. If 
hyperons exist in the interior of NSs (i.e., the nucleonic matter is 
sufficiently stiff), it may be one of the reasons high frequencies 
are not observed. 
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